A fast pace-of-life is traded off against a high thermal performance

The integration of life-history, behavioural and physiological traits into a ‘pace-of-life syndrome’ is a powerful concept in understanding trait variation in nature. Yet, mechanisms maintaining variation in ‘pace-of-life’ are not well understood. We tested whether decreased thermal performance is an energetic cost of a faster pace-of-life. We characterized the pace-of-life of larvae of the damselfly Ischnura elegans from high-latitude and low-latitude regions when reared at 20°C or 24°C in a common-garden experiment, and estimated thermal performance curves for a set of behavioural, physiological and performance traits. Our results confirm a faster pace-of-life (i.e. faster growth and metabolic rate, more active and bold behaviour) in the low-latitude and in warm-reared larvae, and reveal increased maximum performance, Rmax, but not thermal optimum Topt, in low-latitude larvae. Besides a clear pace-of-life syndrome integration at the individual level, larvae also aligned along a ‘cold–hot’ axis. Importantly, a faster pace-of-life correlated negatively with a high thermal performance (i.e. higher Topt for swimming speed, metabolic rate, activity and boldness), which was consistent across latitudes and rearing temperatures. This trade-off, potentially driven by the energetically costly maintenance of a fast pace-of-life, may be an alternative mechanism contributing to the maintenance of variation in pace-of-life within populations.


Introduction
Many organisms exhibit correlated suites of life-history, behavioural and physiological traits forming a pace-of-life syndrome [1,2]. Individuals with a fast paceof-life typically grow and develop faster, show more active, bold and aggressive personalities, have a higher metabolic rate, reproduce earlier, and live shorter [1][2][3]. While proven to be a very useful concept in understanding trait (co)variation in nature, the maintenance of pace-of-life types within populations is still debated [4]. Well-studied trade-offs such as growth versus survival and current versus future reproduction have been suggested as underlying mechanisms of these adaptive trait integration patterns [4,5], hence allowing for the co-existence of individuals with contrasting pace-of-life within a population. Only very recently have researchers started to investigate alternative mechanisms that may explain intra-population variation in the pace-of-life syndrome [6]. Energy-mediated trade-off patterns with tolerance to environmental stressors, such as increasing temperatures [7], are promising candidates in this context, given stress coping mechanisms are energy-demanding [8].
An emerging insight is that organisms can be aligned across a 'cold-hot' axis: individuals situated on the 'hot' end of this axis are said to have a higher thermal performance, i.e. prefer higher environmental temperatures and show a higher maximum performance at higher temperatures, compared to 'cold' individuals (e.g. [9,10]). This axis can be described by reconstructing individual thermal performance curves, describing the relationship between a trait and a temperature gradient [11]. Typically, thermal performance curves have an accelerating rising part (activation energy, E a ) until reaching an optimum temperature (T opt ) at which maximum performance (R max ) is achieved, followed by a fast decelerating part until the critical maximum temperature where performance is zero [11,12]. Individual variation in thermal performance curves reflects differences in the ability to cope with thermal regimes [11,13]. Specifically, the frequently documented individual-level positive correlation between T opt and R max (e.g. [14]) reflects increased performance at higher temperatures, the so-called 'hotter-is-better' hypothesis [12]. Trade-off patterns between a single thermal performance curve parameter and the trait values of a single pace-of-life trait have occasionally been shown. Specifically, a trade-off at the individual level has been detected between a high T opt for locomotor performance and a fast growth rate [15,16]. Yet, such trade-off patterns have not been considered between integrated sets of (i) multiple thermal performance curve parameters of multiple pace-of-life-related traits (i.e. 'cold-hot' axis) and (ii) their trait values (i.e. pace-of-life syndrome). This gap in the literature is important because organismal traits form integrated sets [17], yet their thermal performance curves can be strongly different [11,18]. Furthermore, because environmental conditions can shape trade-off patterns [19], the relationship between pace-of-life traits [7,20] and between thermal performance curve parameters of the same trait [14,21] may depend on the evolutionary background of populations and on experimental treatments (but see: [22]). Such conditional trait integration patterns may also apply to the relationship between the pace-of-life syndrome and thermal performance, yet this has never been tested.
We here tested for alignment at the individual level across a slow-fast (pace-of-life) axis and along a cold-hot (thermal performance) axis, to then evaluate a trade-off relationship between the pace-of-life and thermal performance. We additionally studied whether these among-individual trait integration and trade-off patterns were robust across two latitudes and two rearing temperatures. For this, we conducted a full-factorial common-garden experiment with larvae of the damselfly Ischnura elegans originating from high-latitude (Denmark and Southern Sweden) and low-latitude (Southern France) populations being reared at 20°C or 24°C. We characterized for each larva the pace-of-life by measuring a set of behavioural, life-history and physiological traits, and the thermal performance by estimating thermal performance curves and associated parameters (E a , T opt and R max ) for the same set of behavioural and physiological traits, and for a locomotor performance trait. We explicitly tested at the individual level for pace-of-life and thermal performance patterns, and for a trade-off between these using a structural equation modelling (SEM) framework.
We had predictions both at the level of the latitudes and temperature treatments, and at the individual level. First, based on the higher number of generations per year in the low-versus high-latitude populations (three to four generations per year versus one generation per 2 years: [23]), we expected a faster pace-of-life in the more time-constrained low-latitude populations, as documented in the study species [22,24]. We also expected a faster pace-of-life in larvae when reared at the warmer versus colder temperature, based on the metabolic theory of ecology [25] (e.g. for the water flea Daphnia: [7]). Regarding the thermal performance curves, we expected an overall higher performance, including a higher R max , in low-latitude larvae to cope with the stronger time constraints [26], but not necessarily a higher T opt, given the difficulty of imposing radical changes to the structural and physiological characteristics of evolutionarily conserved enzymes [12]. Nevertheless, a higher T opt in populations originating from warmer latitudes has also been reported and is interpreted as thermal adaptation to local thermal conditions [27,28]. Second, we expected at the individual level the existence of a pace-of-life syndrome as previously shown in the study species [22,29]. Based on intensive work on lizards [9,10], we expected positive covariations between thermal performance curve parameters ( particularly T opt and R max ) of a given trait, as well as across different traits, to form a 'coldhot' axis (e.g. 'hot' individuals have higher T opt and R max for all traits). Finally, given the associated costs of sustaining a fast pace-of-life [1] and a hot thermal performance [13], we expected a negative correlation between the pace-of-life and thermal performance. Hence, we predicted fast-paced larvae to be situated more on the cold end of the cold--hot axis, and vice versa for slow-paced larvae. On the other hand, in case a 'hot' thermal performance or a fast pace-of-life do not carry significant costs, a positive correlation between the pace-of-life and thermal performance may also be expected, as individuals with a proactive behaviour (more aggressive, bold and active) can potentially compensate for the increased metabolic demands at higher body temperatures [9,10].

Material and methods
We studied three low-latitude (southern France) and three high-latitude (Denmark and southern Sweden) populations at shallow lakes within the European range of I. elegans [30] (coordinates in the electronic supplementary material, table S1). We collected egg clutches of ca 10 females per population. We reared larvae individually in 200 ml plastic vials filled with dechlorinated tap water that were placed in incubators at either 20 or 24°C at a photoperiod of 14 : 10 light : dark. Rearing temperatures reflect mean summer water temperatures of shallow lakes in the study region (southern Sweden ca 20°C, southern France ca 24°C: [22,31]). Larvae were fed ad libitum nauplii of Artemia six times a week.
As a general research strategy, we conducted two consecutive sets of trait measurements with distinct aims (electronic supplementary material, figure S1). First, to characterize the paceof-life of each larva, we estimated larval growth rate, activity and boldness behaviour, and metabolic rate at their rearing temperature, in that order. Next, to characterize the thermal performance of each larva, we reconstructed thermal performance curves for the same behavioural (activity and boldness) and physiological (metabolic rate) traits, as well as for swimming speed, at each of six test temperatures (16,20,24,28,32 and 36°C) after a 1 h acclimatisation period. For thermal performance curves, we ran trials at two temperatures per day and conducted trials every other day. We randomized the order of the six test temperatures, except for the highest two temperatures (32 and 36°C), which were conducted on the final trial day to prevent potential carry-over effects of exposure to extreme high temperatures (following: [32,33]). We measured growth rate as gain in body mass over 4 days, activity as total distance moved in 10 min, boldness as latency time to resume activity following a simulated predator attack, metabolic rate as the rate of oxygen consumption, and swimming speed as the speed of the fastest escape swimming bout. Detailed protocols of trait measurements are provided as electronic supplementary material. Sample sizes per treatment combination were 38 high-latitude larvae at rearing temperature 20°C, 53 high-latitude larvae at 24°C, 55 lowlatitude larvae at 20°C and 43 low-latitude larvae at 24°C. For the thermal performance curve reconstructions, this resulted in 1134 measurements for each of four response variables. To test for effects of rearing temperature and latitude on growth rate, activity, boldness and metabolic rate when measured at their rearing temperature ( prior to the thermal performance curve trials), we ran separate linear mixed-effect models (LMMs) with temperature, latitude and their interaction term as fixed effects, and population (nested within latitude) as a random effect. Full statistical results are provided in the electronic supplementary material, table S2. To test for effects of rearing temperature and latitude on thermal performance curves of activity, boldness, metabolic rate and swimming speed, we ran separate LMMs as described above, with the addition of the linear and quadratic terms of test temperature, as well as their interactions with rearing temperature and latitude. To prevent over-parameterization of these quadratic models, we applied stepwise backwards elimination model selection using the step function in the software R. For each trait, we estimated the repeatability [34] across test temperatures (as in [35]), adjusted for fixed effects rearing temperature and latitude (see electronic supplementary material for details of repeatability models).
To test for treatment effects on the extracted thermal performance curve parameters (R max , T opt and E a , see below), we ran separate LMMs with temperature, latitude and their interaction term as fixed effects, and population (nested within latitude) as a random effect. As the random effect of population was never significant (likelihood ratio tests, p > 0.446 for all models described above), we removed this factor from the models. Detailed results of these models (electronic supplementary material, table S3) as well as estimates (with 95% CIs) of thermal performance curve parameters (electronic supplementary material, table S4) are reported as electronic supplementary material.

(i) Estimation of thermal performance curves and their parameters
Following a recently published pipeline for analysing thermal performance curves [36], we fitted a set of candidate thermal performance curve models separately for activity, boldness, metabolic rate and swimming speed, and this for each of the four combinations of rearing temperature and latitude. Based on AIC comparisons, the quadratic thermal performance curve model was found to be the best overall fitting model (details provided in electronic supplementary material). Next, we fitted the quadratic function for each larva separately and extracted three thermal performance curve parameters per trait that are known to have biological and ecological relevance [12]: maximum performance (R max ), thermal optimum (T opt , the temperature at which maximum performance is achieved) and activation energy (E a , the slope of the increasing linear part of thermal performance curves).

(ii) Structural equation modelling
To test the hypotheses that (i) pace-of-life traits integrate to form a slow-fast axis, (ii) thermal performance curve parameters (within and across traits) integrate to form a cold-hot axis and (iii) the slow-fast and cold-hot axes are correlated, we used structural equation modelling (SEM) [37]. SEM is often used in 'syndrome' studies (e.g. [38]), as it allows to construct latent variables that reflect integrated traits sets (here the pace-of-life syndrome and cold-hot axis), and to test for relationships between these latent variables. For this, we constructed SEMs with two latent variables that were hypothesized to separately represent the pace-of-life syndrome and the thermal performance of the larvae (figure 1). The pace-of-life syndrome latent variable comprised four indicator (observed) variables (growth rate, activity, boldness and metabolic rate) that were measured at the rearing temperature (20°C or 24°C) of a larva, prior to the thermal performance curve trials. To refine the models, we excluded several thermal performance curve parameters that were shown to have relatively low predictive power (see electronic supplementary material) and included (residual) correlations between the thermal performance curve parameters and between the pace-of-life traits, as recommended by modification indices (i.e. indices that estimate the amount by which the model χ 2 would be reduced if a single parameter restriction were to be removed from the model).
As a first step, to test the null hypothesis that assumes no trait integration of pace-of-life traits (slow-fast axis) or thermal performance curve parameters (cold-hot axis), we constructed a model with no latent variables. This model was rejected, as indicated by the significant deviation between observed and modelimplied covariances (χ 2 = 135.8, d.f. = 48, p < 0.001), hence was disregarded. Next, to test whether the pace-of-life syndrome and thermal performance were linked, we constructed two models that use the latent structures as described before, whereby one model included a correlation between the two latent variables ('coupled pace-of-life syndrome-thermal performance' model, figure 1a), while the other model did not include this correlation ('independent pace-of-life syndrome and thermal performance' model; figure 1b). As a third alternative hypothesis, we also constructed a non-hierarchical, 'overarching thermal pace-of-life syndrome' model (figure 1c), whereby thermal performance curve parameters and pace-of-life traits were tested for their integration into a single latent variable (as in: [39]). To determine the best-fitting model, we compared the AIC scores of the three royalsocietypublishing.org/journal/rspb Proc. R. Soc. B 289: 20212414 models and further confirmed the fit of the best model (i.e. model with lowest AIC score with ΔAIC > 2) by assessing the fit indices typically used in SEM. Following index values are indicative of models that fit the data well [40]: p > 0.05 for a χ 2 -test, a rootmean-square error of approximation (RMSEA) ≤ 0.05, a standard root-mean-square residual (SRMR) ≤ 0.08 and a comparative fit index (CFI) ≥ 0.95.
Finally, to test whether the integration of pace-of-life syndrome and thermal performance traits depended on the rearing temperature and latitude of origin (or their interaction), we conducted a multi-group SEM using the best-fitting model among the three SEMs (described above) (as in: [20,22]). We used rearing temperature (two levels), latitude (two levels) or the composite variable of rearing temperature × latitude (four levels) as a grouping variable in separate models, and AIC-compared models where the factor loadings on latent variables and correlation between the two latent variables were either the same (constrained model), or allowed to vary between the groups, i.e. treatments (free model). A better fit of the free model would indicate treatment-dependent trait integration patterns.
All analyses were conducted in R v. 4.0.5 [41], using the packages lme4 for LMMs [42], emmeans for post hoc tests [43] and lavaan for SEM [44]. We used the packages rTPC and nls.multstart [36] to fit thermal performance curve models and to calculate thermal performance curve parameters. To conform with model assumptions, we square root-transformed activity and log-transformed boldness.

Results
(a) Effects of latitude and rearing temperature on trait means Growth rates were higher for low-latitude compared to highlatitude larvae (figure 2a; electronic supplementary material,   table S2). Rearing temperature and its interaction with latitude did not affect growth rates. Activity was higher for the low-latitude larvae, and for larvae developed at the higher rearing temperature of 24°C (figure 2b). For boldness behaviour, the significant latitude-by-rearing temperature effect indicated that the higher rearing temperature resulted in increased boldness in the high-latitude (contrast test: p < 0.001) but not in the low-latitude larvae ( p = 0.810) (figure 2c). Metabolic rates were higher for larvae that developed at 24°C compared to 20°C, whereas low-latitude larvae had higher metabolic rates than high-latitude larvae when developed at 20°C ( p = 0.005), but not at 24°C ( p = 0.848) ( p = 0.055; latitude × rearing temperature, figure 2d).
(b) Effects of latitude and rearing temperature on thermal performance curves  3a,b). This was partly mirrored in the slightly higher activation energy (E a ) for activity in low-compared to high-latitude larvae (trend), and in larvae reared at 24°C compared to 20°C (electronic supplementary material,  royalsocietypublishing.org/journal/rspb Proc. R. Soc. B 289: 20212414 except for the T opt being slightly higher in high-latitude larvae (electronic supplementary material, tables S3 and S4). Metabolic rates were overall higher for low-compared to high-latitude larvae, and slightly higher for larvae developed at 20°C compared to 24°C across the test temperatures (table 1 and figure 3e,f ). Higher R max values for low-latitude larvae (electronic supplementary material, table S4) are in line with this pattern, but no significant effect of rearing temperature on R max was detected (electronic supplementary material, table S3). Both the linear increase and the quadratic shape of the thermal performance curve for metabolic rate was affected by the latitude (test temperature × latitude and test temperature 2 × latitude, table 1 and figure 3e,f ). The linear increasing part of the curve was steeper in lowlatitude larvae (estimate ± s.e. = 0.107 ± 0.008) compared to high-latitude larvae (0.075 ± 0.008). This pattern was partially confirmed by the finding of a slightly higher E a (trend) in the low-latitude larvae (electronic supplementary material, tables S3 and S4). Moreover, the thermal performance curve was more concave-shaped, i.e. the quadratic temperature estimate was stronger, for low-latitude larvae (estimate ± s.e. = −0.0018 ± 0.0001) compared to high-latitude larvae (−0.0013 ± 0.0001).
Swimming speeds were overall higher for the lowlatitude larvae across the test temperatures (table 1 and figure 3g,h). Mean swimming speed values, as well as the magnitude of the linear increase and shape of the thermal performance curve, did not depend on any of the treatments (table 1 and figure 3g,h), although the T opt was slightly higher in high-latitude larvae (electronic supplementary material, tables S3 and S4). Using mass-corrected swimming speeds yielded virtually identical results (not shown).

(c) Trait covariation patterns at the individual level
Comparisons of candidate SEMs revealed the 'coupled paceof-life syndrome-thermal performance model' as the best fit for the data (electronic supplementary material, table S7). Note that the other two candidate models were rejected ( p ≤ 0.05, electronic supplementary material, table S7). Rejection of the 'overarching thermal pace-of-life syndrome' model (see figure 1c) indicated that the set of pace-of-life traits and the set of thermal performance traits did not integrate into a single latent variable. This best model (figure 4) indicated (i) the presence of a pace-of-life syndrome latent variable, comprising of growth rate, activity, boldness and metabolic rate (measured at the rearing temperature), (ii) the presence of a thermal performance latent variable, comprising of R max (for swimming speed and boldness), T opt (for swimming speed, metabolic rate, activity and boldness) and E a (for swimming speed and metabolic rate) and (iii) an integration of the two latent variables. The pace-of-life syndrome latent variable indicated that more active larvae were also bolder and had higher metabolic rates and faster growth rates (trend). The two behaviours were the strongest contributors (i.e. had the highest loadings) to the pace-oflife syndrome. As for the thermal performance latent variable, the positive loadings of T opt for all traits indicated that larvae with higher T opt values for one trait also had higher T opt values for the other measured traits. Moreover, larvae with higher T opt values also had a higher R max and E a for swimming speed but lower values of R max for boldness and of E a for metabolic rate. The two strongest contributors to thermal performance were R max and T opt for swimming speed (both with positive loadings). Notably, the two latent variables were negatively correlated, indicating larvae with a faster pace-of-life had lower values for thermal performance (e.g. lower T opt values).
A multi-group SEM approach using AIC scores revealed that the 'constrained' models, i.e. models in which loadings and correlations were fixed across treatments, to have the best fit (see electronic supplementary material, table S8).

Discussion
Next to the findings of a trait integration of life-history, physiology and behavioural traits into a pace-of-life syndrome (indicating a 'slow-fast' axis), and a trait integration of thermal performance traits (indicating a 'coldhot' axis), a major finding of our study was the presence of a trade-off pattern where animals with a faster pace-of-life Table 1. Results of linear models testing for treatment effects on thermal performance curves. t_temp = test temperature (t_temp 2 = quadratic term of test temperature), lat = latitude, r_temp = rearing temperature. Significant effects ( p < 0.05) are indicated in italics. showed a lower thermal performance. To provide background, we first discuss latitudinal and thermal patterns in trait means and thermal performance curves, and then discuss this trade-off.

(a) Pace-of-life across latitudes and rearing temperatures
As expected, low-latitude larvae showed a faster pace-of-life (faster growth and metabolic rate, and higher activity) compared to high-latitude larvae ( [1,45], for the study species: [22]). The higher number of generations produced per year in the low-latitude populations is assumed to be the driver of their faster pace-of-life. Producing multiple generations in a given year requires a fast life-history (i.e. growth rate) which often requires a fast metabolic machinery (i.e. metabolic rate) and a behavioural profile that facilitates increased energy acquisition (i.e. activity). We also confirmed previous findings showing an accelerated pace-of-life (for metabolic rate and activity, and for boldness in high-latitude populations) when reared at a higher temperature (e.g. in the water flea Daphnia magna: [7]). This thermally induced faster pace-of-life was more prominent in high-latitude compared to low-latitude larvae (for metabolic rate and boldness), suggesting a potential compensatory mechanism whereby high-latitude populations more efficiently capitalize on higher temperatures, which they encounter less frequently (see [46] for a similar reasoning in the cooler high-elevation populations of a grasshopper).

(b) Thermal performance curves across latitudes and rearing temperatures
Previous work documenting increased R max in low-latitude populations (e.g. the common woodlouse Porcellio laevis: [27]; the porcelain crab Petrolisthes violaceus: [28]) appears in line with our findings. Yet, these studies also reported higher T opt in lowlatitude populations, in line with the 'hotter-is-better' hypothesis that predicts thermal adaptation to facilitate increased maximum performance at higher temperatures in warm-adapted populations [12,47]. Instead, we observed a higher R max in the absence of a change in T opt in the more time-constrained low-latitude populations. This resembles 'countergradient variation' [26], whereby populations experiencing stronger time constraints show higher (maximum) growth and development rates to compensate for a short growing season. Countergradient variation has, for example, been documented in fish along a latitudinal gradient [48], and in damselflies along an urbanization gradient [49]. The higher maximum performance in the more time-constrained low-latitude populations for activity, boldness and metabolic rate (traits that are tightly linked to an accelerated life history) is in line with this pattern. The mostly similar T opt values across the latitudes indicate no evidence for a thermal adaptation scenario [12]. While higher T opt in the warmer lower latitude populations has been shown in other ectotherms (e.g. [27,28]), the absence of differentiation in T opt across geographically separated populations is not uncommon (e.g. [50,51]) (see also [52] for a review of geographic variation in insect thermal performance curves). Radical changes in the thermokinetics of enzymes, e.g. an increase in temperatures at which an enzyme can function, may require fundamental changes in the genetic architecture, such as mutations in structural genes [12]. Thermokinetics of enzymes are therefore assumed to be evolutionarily conserved systems, which may partially explain findings of relatively stable T opt across latitudes.
(c) Pace-of-life syndrome, cold-hot axis, and their link In accordance with the pace-of-life syndrome hypothesis [1,2], life-history, physiology and behavioural traits showed a positive covariation at the individual level, forming a 'slow-fast' axis. Specifically, larvae at the fast end of the axis had a higher  growth rate, were more active and bolder, and showed a faster metabolic rate. Behaviours that have direct links to energy acquisition, e.g. activity and boldness (for damselflies: [53]), are expected to positively covary with metabolic rate and life history [54,55]. Note that the relatively low trait repeatabilities (all less than 0.25) calculated across the six test temperatures are likely underestimations of the more traditional, across-time repeatabilities, calculated at different time points but at the same temperature. We indeed have observed substantially higher across-time repeatabilities for the study species and for a closely related damselfly species [29,56]. The previous finding of paceof-life syndrome with a similar trait integration pattern [22] further suggests a consistent pace-of-life syndrome in the study system. Next to the slow-fast axis, the tight integration of thermal performance traits at the individual level confirmed the existence of the recently proposed 'cold-hot' axis, that had only been studied in lizards [9,10]. By including thermal performance curve parameters of activity, boldness and metabolic rate, we extended this cold-hot axis, which in previous studies was based on locomotor performance and thermal preference only. Animals at the hot end of the axis showed a higher thermal performance as evidenced by the consistently higher T opt values for swimming speed, metabolic rate, activity and boldness, and also a higher R max for swimming speed, the trait that contributed most to the cold-hot axis. The positive covariation of T opt across traits (see also [11,18]) occurred despite the different thermal performance curves, hence also T opt among traits. The positive association at the individual level between R max and T opt for swimming performance fits the 'hotter-is-better' hypothesis [12] and has been reported previously for locomotor performance (in delicate skinks: [9,10]; in the European common frog: [14]). This pattern was not observed for other traits (metabolic rate and activity) and was even reversed for boldness, suggesting that the thermodynamic constraints assumed to underlie the hotter-is-better hypothesis can be trait dependent.
A key novelty was a coupling of the 'slow-fast' axis and the 'cold--hot' axis, specifically a trade-off pattern where animals with a faster pace-of-life showed a lower thermal performance. Indeed, animals with a faster pace-of-life were situated on the 'colder' side along the cold-hot continuum, i.e. had a lower T opt for all four traits, and a lower R max for swimming speed. The two studies that tested a relationship between (non-behavioural) pace-of-life traits and thermal performance traits only considered locomotor traits, and also reported a negative correlation between growth rate and T opt (tadpoles: [15]; lizards: [16]). Sustaining a faster pace-of-life is assumed to be costly [1] and therefore may require allocation of energy away from costly processes related to thermal performance [13]. The biochemical machinery allowing maximum performance at higher temperatures may indeed be traded off with a fast pace-of-life. For example, an increase in maximum thermal tolerance, a thermal performance parameter often shown to positively covary and coevolve with T opt [57], has been shown to trade-off with growth rate (e.g. in the rainbow trout: [58]). Another potential mechanism underlying a high T opt , namely the costly production of heat shock proteins to prevent protein denaturation at high temperatures [59], has also been shown to trade-off with growth rate [60] (for the study species: [61]). Whatever the exact underpinnings of the negative coupling between the pace-oflife syndrome and the cold-hot axis, this trade-off pattern may be an ignored, but potentially important mechanism for the co-existence of individuals with contrasting pace-of-life within a population. The mechanisms maintaining variation of the pace-of-life within a population are still highly debated [1,4] and are typically based on trade-offs between pace-of-life traits; i.e. growth versus survival [5], or current versus future reproduction [4]. Our data suggest that another reason why fast-paced individuals may not always be favoured is that they have a lower thermal performance.
Intriguingly, the two other studies that explicitly looked at the link between behavioural pace-of-life traits and thermal performance, both on the terrestrial skink Lampropholis delicata, found instead a positive covariation whereby animals showing higher exploration, activity and boldness had also a higher T opt and R max for sprint speed, and a higher thermal preference [9,10]. These authors therefore suggested to extend the traditional pace-of-life syndrome to include thermal performance. The absence of a trade-off between fast pace-oflife and thermal performance may point to reduced costs of having a 'hot' thermal type in this system. Potentially, because lizards are commonly exposed to strongly fluctuating, hence also very high temperatures, they may be energetically adapted to develop 'hotter' thermal types with reduced energy allocation. It is important to note that such a positive association between behaviour and thermal performance may indeed be expected for terrestrial, heliothermic animals that rely heavily on behavioural thermoregulation such as basking behaviour, as lizards with a 'hot' type that also have a more proactive behavioural type (i.e. more active, exploratory and bold) can generate more basking opportunities. Given that aquatic insects such as damselfly larvae are better protected from drastic temperature fluctuations (both in time and space) due to the buffering effect of water, a positive relationship between pace-of-life and thermal performance driven by thermoregulatory behaviour is much less expected. A recent study indeed reported no link between preferred temperature and behavioural type, i.e. bold versus shy, in a freshwater fish [62] (but see [63]).
Another key finding was that the covariation patterns, both among pace-of-life traits and among thermal performance traits (i.e. thermal performance curve parameters) as well as the trade-off pattern between pace-of-life and thermal performance, were stable across latitudes and rearing temperatures. Consistency in the structure of pace-of-life syndromes (see also: [22]) and thermal performance traits may impose restrictions on independent (evolutionary) responses of individual traits to warming, assuming the syndrome structure is underpinned by genetic correlations [64]. Yet, note that this may not necessarily constrain adaptive evolution of the integrated trait sets (i.e. the pace-of-life syndrome and thermal performance), given that adaptive trait integration may potentially guide adaptive evolution of sets of traits along the slow-fast continuum. Notably, also the negative relationship between pace-of-life and thermal performance was robust against potential influences of rearing temperatures and latitude of origin, suggesting that the trade-off between a fast pace-oflife and thermal performance is general in the study species.
In summary, our study extends the recent finding of an integration between behavioural traits and thermal traits for locomotor performance [9,10] by including life history and physiology; traits typically studied under the pace-of-life syndrome hypothesis. We thereby identified a novel trade-royalsocietypublishing.org/journal/rspb Proc. R. Soc. B 289: 20212414 off pattern, revealing fast-paced individuals to have lower optimum temperatures, potentially due to the energetically costly maintenance of a fast pace-of-life and/or of high optimal temperatures. We propose this trade-off to be an alternative mechanism that may contribute to the maintenance of variation in pace-of-life within populations; a topic that is still unresolved [1,4].